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Abstract. Here we investigate to what extent unresolved 7 -ray pulsars contribute 
to the galactic diffuse emission, and further whether unresolved 7 -ray pulsars can be 
made responsible for the excess of diffuse galactic emission above 1 GeV which has 
been observed by EGRET. Our analysis is based only on the properties of the six 
pulsars which have been identified in the EGRET data, and is independent of choice 
of a pulsar emission model. 

We find that pulsars contribute very little to the diffuse emission at lower energies, 
whereas above 1 GeV they can account for 25% of the observed intensity in selected 
regions for a reasonable number of directly observable 7 -ray pulsars (~12). While the 
excess above 1 GeV 7 -ray energy is observed at least up to six or eight degrees off the 
plane, the pulsar contribution would be negligible there. Thus pulsars do significantly 
contribute to the diffuse galactic 7 -ray emission above 1 GeV, but they can not be made 
responsible for all the discrepancy between observed intensity and model predictions 
in this energy range. 


A recent analysis of the diffuse galactic 7 -ray emission in the energy range of 30 
MeV to 30 GeV [1] indicated that the spatial structure and total intensity of the 
emission observed by EGRET can be well understood as the result of interactions 
between cosmic rays and the interstellar medium in addition to a isotropic extra- 
galactic background. However, at energies above 1 GeV the models predict only 
roughly 60% of the observed intensity. 

The attempt of this paper is twofold. We want to provide constraints on the 
general contribution of the most likely input from discrete sources - pulsars - to the 
diffuse galactic 7 -ray emission. We also want to find out whether or not unresolved 
pulsars can account for the observed excess at high 7 -ray energies. Instead of using 
models as was done in previous studies [2-4] we will base our analysis solely on 
the properties of the six pulsars observed by EGRET. A detailed description of the 
method can be found elsewhere [5]. 

EGRET has up to now identified six pulsars by their light curves. The differential 


7 -ray photon spectra of these sources have been derived from pulsed analysis [ 6 ]. 
We will concentrate on nine energy bands spanning 50 MeV to 10 GeV, for which 
we have the observed photon flux Sj(E k ) of pulsar i in energy band E k . The basis 
of our modelling is the assumption that pulsars do not behave arbitrarily in their 
relation between the distance-normalized y-ray intensity and the age, but that they 
follow a trend, a correlation between intensity and age. The best evidence comes 
from the data themselves: the y 2 sums we obtain indicate that a correlation is 
an appropriate description of the data. The limited number of degrees of freedom 
forces us to use the most simple correlation model, a power-law: 


S m {t,E k ) = 10 yk (—— - \ kpc 2 /cm 2 /sec/MeV 
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This relation can be fitted to the data on the basis of weighted least squares for 
all energy bins. The distance Di and its uncertainty can be generally derived on 
the basis of the pulsar dispersion and rotation measure [7] and only for extremely 
nearby objects like Geminga by parallax measurements [ 8 ]. The age of pulsars can 
be estimated from the ratio of period and period derivative. The uncertainty of 
this measure of age is large, especially since pulsars often do not slow down solely 
by dipole radiation. We will use a factor of two for the age uncertainty St t , except 
for the Crab for which the true age is taken with a nominal uncertainty of 10 %. 

Three sources of uncertainty have to be considered in the fit: uncertainties in 
intensity, in distance, and in age. While the intensity error can be assumed to 
follow a Gaussian probability distribution, both the errors in age and distance enter 
with some power, so that their effective probability distributions are definitely not 
Gaussian. Furthermore, the way how age and distance estimates are derived lets 
us think that a Gaussian probability function is not a fair description of the actual 
error distribution of age and distance, respectively. We can account for these effects 
and still use the y 2 method, if we assume that the error distribution for age and 
distance are Gaussians in the logarithm, i.e. the uncertainty is multiplicative rather 
than additive. Only in this case does any power of these parameters obey the same 
corresponding error distribution. 

To account for the different uncertainty distributions in intensity on one side 
and age and distance on the other side, we have separated each argument \i °f 

the y 2 -summation into three components y 2 = (y % 2 + Xl‘i + Xi,s) > one f° r each 
source of uncertainty, where the components can be calculated in the system in 
which the uncertainty distribution may be taken as Gaussian. The best fit values 
and uncertainties of the parameters y k and b k of Eq.l will serve as input for the 
calculation of the pulsar contribution to the diffuse galactic y-ray emission. 

For each source we can define a critical distance, up to which the source can be 
detected directly, and beyond which it would contribute to the diffuse emission. 
We will assume that the Galaxy is a simple disk-like entity of radius ry and half¬ 
thickness Zh, in which the line-of-sight through the disk determines the sensitivity 
threshold F c (l,b). The numbers ry, Zh have been chosen such that on average the 
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FIGURE 1 . Here we show the fraction of the total diffuse intensity in direction of the Galactic 
Center (left) and galactic anticenter (right) which can be attributed to unresolved pulsars. The 
parameters to the model are given in the text. The error bars are derived by propagation of the 
parameter uncertainty in the pulsar model fit. The intensity due to pulsars scales almost linearly 
with the number of objects which are supposed to be seen as point sources, in this case 11.9 
pulsars. 

threshold corresponds to the flux of a 5 cr source, respectively 4 a for |6| >10°, in the 
summed data of Phases 1-4. The threshold in direction of the galactic poles would 
be F c (0,90) = 7 • 10 ” 8 ph./cm 2 /sec, whereas in direction of the galactic center it 
would be F c (0,0) = 4.8 • 10”' ph./cm 2 /sec. Then the critical distance can to be 
calculated as 
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where 8E k is the width of the energy bin E k . The spatial distribution of pulsars 
is poorly determined. The distribution of radio pulsars indicates that there is 
a galactocentric gradient [7]. We may thus parametrize the normalized spatial 
distribution of pulsars in galactocentric cylinder coordinates (r,z) as 

p(r,z) = 0.0435 @(| z c — z\) (^ c r ^) _1 cosh ” 1 (— 'j (3) 
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where the cosh-term accounts for the galactocentric gradient and 0 is a step func¬ 
tion. With Tp 1 as the birth rate of y-ray pulsars and e as the fraction which radiates 
in our direction we can determine the number of directly detectable pulsars by in¬ 
tegration over the line-of-sight, solid angle, and age 
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FIGURE 2. The latitude distribution of the diffuse emission above 1 GeV as observed by 
EGRET is much wider than that of the pulsar contribution, here shown as solid line for 12 
directly observable pulsars. Both the data and the model have been averaged over longitude. 
Here we show the original observed intensity distribution compared to the psf-convolved model 
prediction. The statistical error of the observed emission is at the percent level and thus negligible. 


as well as the diffuse emission of unresolved pulsars 

f rtmax r POO 

I dif {E k ) = —- dt I dil' dx p(r,z) S m (t,E k ) (5) 

Tp \ L J 0 J ^2 J Xmax 

where x is the distance coordinate along the line-of-sight. Please note that any 
change in the total number of pulsars has similar impact on the number of directly 
observable pulsars as on the y-ray intensity of unresolved objects. Thus the direct 
detections provide a strong constraint on the pulsar contribution to the diffuse 
galactic y-ray emission and limit our choice of e/r p . With Eq.5 we can now calculate 
the spectrum of unresolved pulsars and compare it to the spectrum of the total 
observed diffuse emission. This is shown for the Galactic Center direction and for 
the anticenter in Fig.l, where we plot the fraction of the total observed emission 
which is due to pulsars. Here we used for the spatial distribution of pulsars a 
radial scale length r c =3.5 kpc and a half-thickness of z c =0.15 kpc. The result 
changes little when we choose an age-dependent thickness which would account for 
older pulsars being more distant from the plane. The maximum age t max , we have 
used 2 ■ 10 6 years, has only little influence on the result since it affects the number 
of directly observable objects N det nearly as much as the predicted emission of 
unresolved pulsars I,j x / • 

The parameters used here (r -1 = 0.01 years -1 and e = 0.15) imply that around 
12 pulsars should be detectable by EGRET. Since six are already identified this 




would mean that another six unidentified sources are actually pulsars. The small 
fraction of unresolved EGRET sources which shows pulsar-like or Geminga-like 7 - 
ray spectra [9] argues strongly against the bulk of unidentified sources being pulsars. 
Also a substantial fraction of the unidentified EGRET sources in the galactic plane 
appears to be variable which makes the identification of all unidentified sources 
with pulsars even more unlikely [10]. In total, we think that around 12 directly 
observable pulsars, of which six are already identified, is a reasonable number, and 
that thus the calculated contribution of unresolved pulsars can be taken as serious 
estimate. Considering the integrated intensity above 100 MeV in this example pul¬ 
sars would cause 13% of the observed intensity in the Galactic Center direction and 
about 5% in the anticenter direction. We have compared the latitude distribution 
of the observed diffuse emission above 1 GeV to what our model predicts for the 
pulsar contribution. The result is given in Fig.2 where we show the distribution 
of the diffuse intensity at energies above 1 GeV as averages over longitude. We 
find that pulsars contribute only very close to the plane where the line-of-sights are 
long, but the observed emission does fall off much less rapidly with latitude than 
the emission of unresolved pulsars. Also the observed spectral discrepancy above 
1 GeV seems to extend to latitudes | 6 | > 5°. Thus there must be additional effects 
playing a role for the observed diffuse 7 -ray emission above 1 GeV. 
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